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ABSTRACT 

We show that the gas in growing density perturbations is vulnerable to the influence of winds out- 
flowing from nearby collapsed galaxies that have already formed stars. This suggests that the formation 
of nearby galaxies with masses <; 10 9 M Q is likely to be suppressed, irrespective of the details of galaxy 
formation. An impinging wind may shock heat the gas of a nearby perturbation to above the virial 
temperature, thereby mechanically evaporating the gas, or the baryons may be stripped from the per- 
turbation entirely if they are accelerated to above the escape velocity. We show that baryonic stripping 
is the most effective of these two processes, because shock-heated clouds that are too large to be stripped 
are able to radiatively cool within a sound-crossing time, limiting evaporation. The intergalactic medium 
temperatures and star-formation rates required for outflows to have a significant influence on the forma- 
tion of low-mass galaxies are consistent with current observations, but may soon be examined directly 
via associated distortions in the cosmic microwave background, and with near-infrared observations from 
the Next Generation Space Telescope, which may detect the supernovae from early-forming stars. 

Subject headings: galaxies: formation - intergalactic medium - cosmology: theory 



1. INTRODUCTION 

Outflows from dwarf starbusting galaxies have long been 
the subject of theoretical and observational investigation. 
Since the 1970s, theoretical work has shown that super- 
novae (SNe) and OB winds in these low-mass objects pro- 
duce energetic outflows that shock and enrich the inter- 
galactic medium (IGM) in which they are embedded (Lar- 
son 1974; Dekel & Silk 1986; Vader 1986). This behavior 
has been clearly identified in studies of both local star- 
bursting galaxies (Axon & Taylor 1978; Heckman 1997; 
Martin 1998) and their environments (dclla Ceca et al. 
1996; Bomans, Chu, & Hopp 1997) as well as with spec- 
troscopy of high-z galaxies (Franx et al. 1997; Pettini et 
al. 1998; Frye, Broadhurst & Sprinrad 1999). 

Whether these outflows lead to a catastrophic loss of 
the interstellar gas in these objects is much more uncer- 
tain, with recent simulations and observations suggest- 
ing that dwarfs are very inefficient in removing gas from 
their cores (Mac Low & Ferrara 1999; Murakami & Babul 
1999). However, the issue of whether dwarfs retain a size- 
able fraction of their gas is to some degree decoupled from 
the formation of outflows. In galaxies in the mass range 
1O 7 M ^ M <; 1O 9 M , Mac Low & Ferrara (1999) have 
shown that a "blowout" occurs in which the superbub- 
bles produced by multiple SNII explosions punch out of 
the galaxy, shocking the surrounding IGM while failing to 
excavate the interstellar medium of the galaxy as a whole. 

In the generally investigated hierarchical models of 
structure formation, low-mass galaxies form in large num- 
bers at early times (e.g. White & Frenk 1991) are highly 
clustered (Kaiser 1984). The existence of a large number 
of small and clustered galaxies at early times is also fa- 
vored observationally by the steep number counts and low 
luminosities of faint galaxies (Broadhurst, Ellis & Glaze- 



brook 1992), the sizes of faint galaxies in the Hubble Deep 
Field (HDF) (Bouwens, Broadhurst, & Silk 1999a,b), and 
the clustering properties of Lyman-break galaxies (Adel- 
berger et al. 1998). Locally however, the number density 
of low-mass objects is far less than predicted theoretically 
(Ferguson & Binggeli 1994; Klypin et al. 1999; Moore et al. 
1999), leading to theoretical studies into the disruption of 
dwarves by tidal forces from neighboring objects (Moore, 
Lake, & Katz 1998) and external UV radiation ( Norman 
& Spaans 1997; Corbelli, Galli, & Palla 1997; Ferrara & 
Tolstoy 2000). 

Little attention has been directed however, toward the 
influence of outflows on the formation of neighboring 
galaxies, despite the fact that typical outflow velocities 
exceed the virial velocities of objects in the mass range in 
which the overabundance is most severe. In a companion 
article (Scannapieco & Broadhurst 2000) we describe the 
effects of heating and enrichment by outflows on galaxy 
formation using a Monte Carlo treatment of hierarchical 
structure formation. Here we show using simple scaling 
arguments that suppression of galaxy formation by out- 
flows from dwarf galaxies is important over a large range 
of halo masses, irrespective of the details of cosmological 
simulations. 

The structure of this work is as follows. In §2 we con- 
sider two scenarios for the suppression of dwarf galaxy 
formation by outflows. In §3 we show that the outflow 
models considered fall well within the bounds of current 
observations. Conclusions are given in §4. 

2. SUPPRESSION OF GALAXY FORMATION 



We consider two processes by which outflows from neigh- 
boring dwarves inhibit the formation of a galaxy. In 
the "mechanical evaporation" scenario, the gas associated 
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with an overdense region is heated by a shock to above its 
virial temperature. The thermal pressure of the gas then 
overcomes the dark matter potential and the gas expands 
out of the halo, preventing galaxy formation. In this case, 
the cooling time of the collapsing cloud must be shorter 
than its sound crossing time, otherwise the gas will cool 
before it expands out of the gravitational well and will 
continue to collapse. 

Alternatively, the gas may be stripped from a pertur- 
bation by a shock from a nearby source. In this case, the 
momentum of the shock is sufficient to carry with it the 
gas associated with the neighbor, thus emptying the halo 
of its baryons and preventing a galaxy from forming. In 
the following we evaluate the importance of these two ef- 
fects in turn. 

2.1. Mechanical Evaporation 

The first mechanism we consider is the shock heating of 
the halo gas to a temperature sufficient to cause the gas 
to evaporate into the IGM. This will happen for cases in 
which the virial temperature of the object T v is less than 
the postshock temperature T s , where T s can be written for 
an adiabatic shock as 
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where p is the molecular weight of the gas, which we take 
to be 0.6, m p is the proton mass, k is Boltzmann's con- 
stant, and v s is the velocity of the expanding shock. 

For the typical outflows we will be considering, the 
blast wave reaches the neighboring perturbation of in- 
terest within one cooling time t c « where n e = 

4.5 x 10~ 7 (1 + z) 3 h 2 cm -3 is the average number density 
of electrons, A is the cooling function which we normalize 
to A_22 = A/10~ 22 erg cm 3 s~\ and h is the Hubble con- 
stant in units of 100 km/s/Mpc. Here and below we take 
the baryonic density in units of the critical density to be 
Qb = 0.05. Thus the shock velocity can be approximated 
by a Sedov- Taylor blast wave solution, 



v s = 1600 (eAf/i) 1/2 ^cW 



(2) 



where Af is the number of supernovae driving the bubble 
(each with a total energy 2 x 10 51 erg to take into ac- 
count the contribution from stellar winds), e is the total-to- 
kinetic energy conversion efficiency, <i Cj kpc is the comoving 
distance from the explosion site in kpc xh. The Compton 
drag on the expanding shock as well as the Hubble ex- 
pansion are easily shown to be negligible for our purposes 
in this paper (Ferrara 1998; Scannapieco & Broadhurst 

2000). From the relation T v = 70M 2/3 (1 + z)f^ /3 K, 
where M 6 = M/(W 6 M Q /h) and f2 is the matter density 
in units of the critical density, it follows that collapsed 
objects with total masses lower than M$ < 3.5 x 10 8 (1 + 

z)- 3 / 2 (eAfh) 3 / 2 d~% 2 c n Q 1/2 will be affected by shock heat- 
ing and evaporation. We estimate d Cj kpc as the mean spac- 
ing between objects of a mass scale M, which is appropri- 
ate as the outflowing objects and the forming density per- 
turbations of interest are of roughly the same mass scales: 
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M 6 £ 41(1 + z)-^ 5 (eAfh)^ 5 n 2 /5 . 
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Assuming that one SN occurs for every 100 Mq of 
baryons which form stars (see eg. Gibson 1997) we can 
relate the number of SNe and the mass of the halo as 
Afh = 500Mge s / S1q , where e s / is the initial star forma- 
tion efficiency. Assuming an initial star formation effi- 
ciency of e s f = 0.1 and a kinetic energy conversion effi- 
ciency of e = 0.2 we can estimate eAfh as approximately 
IOMqQq , ano - t nus we ex P ec t d\fh ~ SOOOSIq" 1 for typical 
dwarf-galaxy sized objects of ~ 5 x 10 8 Mq. 

The relation given by Eq. (Q) is shown by the upper 
lines in Fig. 1 for ej\fh = 10 4 . Also shown for reference on 
this plot is the mass below which effects due to photoevap- 
oration are important (M 6 <; 4.4 x 10 3 (1 + z)^ 3 / 2 Q^ 1 ^ 2 ). 
Galaxies with masses below this limit are readily evapo- 
rated by the UV background after the reionization of the 
universe (Barkana & Loeb 1999; Ferrara & Tolstoy 2000) 

2.1.1. Cooling 

To suppress galaxy formation, however, it is not enough 
only to heat the halo gas, because the gas could cool 
rapidly and fall back towards the center of the object. 
Therefore Eq. is a necessary but not sufficient con- 
dition for mechanical evaporation. The second condition 
requires that the gas remains hot for the time required to 
escape the system, or that the cooling time be longer than 
the the sound crossing time, ^4- > (/c s , where I is size 

of the halo, c s = (kT s / fimj,) 1 ' 2 is the sound speed, and 
n e = n e 8 where 8 = p/ po is the overdensity of the region, 
po is the mean density of the universe. This gives 



M 6 <; 2.3(dV/i) 9/11 (l + z)- 12/11 ^ 4/11 /i- 6/n 
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which forms the upper bound of the cross-hatched region 
in Figure 1. Here we see that for both cosmological models, 
the cooling vastly undercuts the minimum mass necessary 
for disruption by heating, dropping below the photoevap- 
oration limit in the ACDM model. While it is relatively 
easy for galactic outflows to shock low-density gas to above 
the virial temperatures of the potential, this excess energy 
is efficiently radiated away for most large clouds. 

2.2. Baryonic Stripping 

The second possibility is that all the gas is striped by 
the impinging shock. We may estimate that this will occur 
when the shock moves through the center of the halo with 
sufficient momentum to accelerate the gas to the escape 
velocity of the halo: fM s v s > M c v e where / = £ 2 /4d 2 
is the solid angle of the shell impinging on the lump, M s 
is the mass of material swept up by the expanding shock, 
M c = (4ir/3)p c £l},S£ 3 is the baryonic mass of the cloud 
whose radius is £, and v e — (GMjl) 1 ! 2 is the escape speed 
from the cloud. Solving for the mass that can be swept up 
as function of redshift and replacing d c ^ pc by the mean 
separation as estimated from Eq. (|3|) gives 



M 6 £ 53<r x (l + z)- 3 ^ 5 (eAfh) 3/5 n 2 /5 . 



(6) 
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which is also plotted in Figure 1. Here we see that while it 
is significantly easier to heat a cloud than to sweep away 
its baryons, the short cooling times of most halos cause 
baryonic stripping to have the greatest impact on the for- 
mation of nearby galaxies. Note that in this plot the den- 
sity is taken to be only twice that of the mean density of 
the universe. Thus it is important that the perturbation 
be in the linear or weakly nonlinear regime for this mecha- 
nism to be effective. As the time between turn-around and 
virialization is relatively short, however, most protogalax- 
ies are likely to be at these low density contrasts when 
shocked by neighboring objects. 

3. OBSERVATIONAL CONSTRAINTS 

3.1. Compton y-parameter 

While widespread baryonic stripping of dwarf galaxies is 
a generic consequence of associating outflows typical of lo- 
cal dwarfs with primordial dwarf galaxies, it is important 
to check that this extrapolation to higher rcdshift is con- 
sistent with observational constraints. An unavoidable re- 
sult of widespread IGM heating is the presence of spectral 
distortions in the cosmic microwave background (CMB). 
The degree of these distortions is given by the Compton y 
parameter which is simply the convolution of the optical 
depth with the electron temperature along the line of sight 
(Zel'dovich & Sunyaev 1969; Sunyaev & Zel'dovich 1972). 
To calculate the mean y we make use of the fact that the 
total optical depth within ionized regions must be within 
the observational limit of r < 0.5 (Griffiths, Barbosa, & 
Liddle 1999). As the total optical depth within outflows 
must be less than the total optical depth overall we can 
estimate this convolution as 

Kffr^), (7) 

14 \J r . m e c 2 J 

where r.; is the initial size of the blast, which we take to be 
l/i -1 comoving kpc, and V4 is the volume defined by the 
radius outside which the temperature of the blast drops be- 
low the ionization temperature of hydrogen w 10 4 , which 
is readily estimated from Eq. (Q) and (0) to be equal to 
15 (eA/Ti) 1 / 3 h- 1 (1 + z)- 1 kpc. 
This gives 

/f 1 4 K 

y&T^-^MlbieAfh^zlAxW- 5 , (8) 
m e c z 

where we take (eJ\fh) = 10 4 . This is within the bounds set 
by the COBE data of y < 1.5 x 1(T 5 (Fixsen et al. 1996), 
but only marginally. This suggests that future CMB ex- 
periments may help to constrain the degree of IGM heating 
by outflows. 

3.2. Point Source Luminosities 

A second consequence of widespread dwarf outflows is 
the presence of objects seen as point sources in deep sur- 
veys. A SN number of Af « 5 x 10 4 implies that the 
local star formation rate is roughly 0.5M Q /yr. The ques- 
tion arises if objects with these luminosities could have es- 
caped detection as point sources in the HDF. Assuming a 
Salpeter initial mass function, Ciardi et al. (2000) give the 
luminosity per unit frequency at the Lyman limit obtained 



from the adopted spectrum of a primordial galaxy at early 
evolutionary times as jo = 4 x 10 20 M* ergs _1 Hz _1 , where 
M* is the stellar mass of the object. Scaling to the stellar 
mass value implied by the assumed efficiency of massive 
star formation, M* = 5 x lO 6 A/ , we find for the luminos- 
ity of the SN parent cluster of stars L u ~ 6 x 10 42 erg/s 
at the Lyman limit. The observed flux is then J-(vq) = 
L„(l + z)/47rd 2 , where v$ = + z) is the observed fre- 
quency. This gives approximately T = 12(1 + z)~ 1 h 2 nJy 
In order to be observed in the HDF optical filter centered 
at 6060A (V606), the object should be located at redshift 
1 + z = 6.6, thus with a flux equal to ss 2h 2 nJy. This value 
corresponds to a AB magnitude Vab — 30.6 which is well 
below the limiting magnitude of the HDF for point sources 
(typically 28 mag, see Mannucci & Ferrara 1999). How- 
ever, these objects may be detectable by the Next Genera- 
tion Space Telescope, which should reach AB magnitudes 
of order 32 in the near infrared, although dust extinction 
may complicate this analysis. 

4. CONCLUSIONS 

Using simple scaling arguments, we have shown that 
outflows from low-mass galaxies play an important role in 
the suppression of the formation of nearby dwarf galaxies. 
While the details of this process depend on the spatial dis- 
tribution of forming galaxies, and can be only be studied 
in the context of cosmological simulations (Scannapieco & 
Broadhurst 2000), the overall implications of this mecha- 
nism can be understood from the general perspective of 
hierarchical structure formation. 

In principle, outflows can suppress the formation of 
nearby galaxies both by shock heating and evaporation 
and by stripping of the baryonic matter from collapsing 
dark mater halos; in practice, the short cooling times for 
most dwarf-scale collapsing objects suggest that the bary- 
onic stripping scenario is almost always dominant. This 
mechanism has the largest impact in forming dwarves in 
the <^ 10 9 M Q range which is sufficiently large to resist 
photoevaporation by UV radiation, but too small to avoid 
being swept up by nearby dwarf outflows. 

It is interesting to note that numerical studies working 
from a different perspective have similarly shown that mo- 
mentum transfer is the most important feedback mecha- 
nism in determining the internal structure of larger galax- 
ies. While SN heating is relatively ineffective at regulat- 
ing star formation in numerical models, accounting for the 
kinetic component of SN feedback has resulted in quasi- 
stationary models of disc galaxies with reasonable star for- 
mation rates (see eg. Katz 1992; Navarro & White 1993; 
Mihos & Hernquist 1994; Springel 2000). 

Various analytical and N-body studies (Kauffmann, 
White, & Guiderdoni 1993; Klypin et al. 1999; Moore et 
al. 1999) have shown that ~ 50 satellites with circular ve- 
locities ~ 20 km/s and masses ~ 10 9 M Q should be found 
within 600 kpc of the Galaxy, while only 11 are observed. 
While the fraction of objects suppressed is dependent on 
the assumed cosmology and spatial distribution of galax- 
ies, the partial suppression of the formation of objects in 
this mass range is a general consequence of baryonic strip- 
ping by outflows. Note that dark matter halos which are 
subject to sweeping are likely to accrete some gas at later 
times and thus this scenario provides a natural mechanism 
for the formation of "dark" Milky Way satellites, which 
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may be associated the abundant High- Velocity Clouds as 
discussed by Blitz et al. (1999). Also, this population of 
halos could be identified with the low mass tail distribu- 
tion of the dark galaxies that reveal their presence through 
gravitational lensing of quasars (Hawkins 1997). 

The existence of an era of widespread IGM heating 
through outflows is consistent with current CMB con- 
straints, but causes a mean spectral distortion in the range 
that will be probed by the next generation of experiments. 
Similarly, the luminosities of the majority of the starburst- 
ing dwarfs that contribute to this process are beyond the 



limiting magnitudes of the Hubble Deep Field, but may be 
observable with of the Next Generation Space Telescope. 
Thus it may be in the near future that we can directly 
observe the era in which dwarf galaxies were first formed 
and the impact of galactic outflows on this process. 

We would like to thank Marc Davis, Joseph Silk, and an 
anonymous referee for helpful comments and discussions. 
This research was supported in part by the National Sci- 
ence Foundation under Grant PHY94-07194. 
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Fig. 1. — Relevant mass scales for suppression of dwarf galaxy formation. The upper lines are the masses below which halos will be heated 
beyond their virial temperatures, although cooling prevents mechanical evaporation from occurring for halos with masses above the cross 
hatched regions. The second highest set of lines, bounding the lightly shaded regions, show the masses below which baryonic stripping is 
effective. Finally the heavily shaded regions show objects that are susceptible to photoevaporation. The upper panel is a flat CDM model and 
the lower panel is a flat ACDM model with Qq = 0.3. In all cases (eNh) = 5000Q \ the overdensity <5 = p/p = 2.0, tt b = 0.05, h = .65, and 
the cooling function A_22 = 1 as described in §2.1.1. Note that photoevaporation affects a larger mass range than mechanical evaporation in 
the ACDM cosmology. 



